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Abstract. For the 2d Euler dynamics of patches, we investigate the conver- 
gence to the singular stationary solution in the presence of a regular strain. 
We prove that the rate of merging can be made double exponential. 



< 

^ ■ 1. Introduction and statement of the results 

> ■ 

In this paper, we study the 2d Euler dynamics of patches: 

= V0 • u, u = \7 ± A~ 1 9, 9{z, 0) = 6 , z= (x, y) £ R 2 , V x = (-d y , d x ) (1) 

where the 9q is the characteristic function of some compact set (a patch). In our 
case, this compact set will be a centrally symmetric pair of two simply connected 
1 domains with rather regular boundaries. 

This problem attracted a lot of attention in both physics and mathematical lit- 
" erature in the last several decades and became a classical one. The existence of 

<n: weak solutions for Euler equation is due to Yudovich |12j ; his result, in particular, 

ensures that the dynamics of patches is well-defined. In [3], Chemin proved that if 
the boundary of the patch is sufficiently regular then it will retain the same regu- 
larity forever; another proof of that fact was given later by Bertozzi and Constantin 
P]. We recommend the wonderful books [2111] for introduction to the subject and 
for simplified proofs. For closely related models (e.g., SQG), there were attempts 
Jv> . recently to prove that a singularity can occur in finite time (see, e.g., [5] and related 

In; UEW- 

Although we have global regularity for the 2d Euler dynamics of contours, very 
little is known about the mechanism leading to "singularity formation" . The latter 
can be described in various ways: growth in time of the curvature, perimeter, etc. 
In this paper, we focus on one particular geometric characteristic, the distance 
between two interacting patches, and show that it can decay as double exponential 
infinitely in time. We also reveal the mechanism that drives this process. 

Let f2~ = { — z,z G fi}, i.e., the image of f2 under the central symmetry. The 
boundary of will be denoted by T. 

In R 2 ~ C, we consider the 2d Euler dynamics given by the initial configuration 
f2(0)Uf2 _ (0). The areas of patches, the distance between them, the value of vorticity 
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- all these quantities are of order one as t = 0. As time evolves, the Euler evolution 
deforms fl(0) U fi - (0) to a new pair fi(i) U Q~(t) with \Q(t)\ = |fl(0)| because 
the flow preserves the area and the central symmetry. These new patches will be 
separated from each other for all times, i.e. ft(t) n f2~(i) = 0, but the question, 
however, is how small the distance between them can get? The main result of the 
paper is 

Theorem 1.1. Let positive 5i be sufficiently small. Then, there is a simply con- 
nected domain Q(0) satisfying dist(fi(0), f2 — (0)) ~ 1, and a time- dependent com- 
pactly supported incompressible odd strain 

S{z,t) = (P(z,t),Q(z,t)) 

such that 

dist(ft(i),0) ~ dist(fi(i),fi _ (f)) < exp(-e ,5lt ) 
where f2(t) U is the Euler dynamics of f2(0) U f2 _ (0) in the presence of the 

strain S(z,t). Moreover, for S(z,t) we have 

\S(z,t)\ 

sup 1 , n < oo (2) 

t>0,2#0 \z\ 

Remark 1. The corresponding equation for 9{z,t) = Xn(t) + Xn-(t) is 

9 = V9 ■ (u + S), u = V 1 - A~% 6(z, 0) = xn(o) + Xa- (o) (3) 

(compare this to {TJ). 

Remark 2. As it will be clear from the proof, these contours will touch each 
other at t = +oo and the touching point is at the origin. In the local coordinates 
around the origin the functions parameterizing the contours converge to ±|x| in 
a self-similar way which will be described in detail. In the following elementary 
lemma, we show that under the S'-strain alone no point can approach the origin 
in the rate faster than exponential as long as assumption @ is made and so it 
is the nonlinear term V6* ■ u in ([3]) that produces the "double exponentially" fast 
singularity formation. 

Lemma 1.1. Let S(z,t) be an odd vector field that satisfies Consider 9 {z,t) — 
Xr(t)( z ) which solves 

6 = V0-S(z,t), 0(z,O)=xt(o)(*) (4) 

where dist(Y(0),0) > and T(0) is a compact set. Then, 

dist(T(t),0) > e~ ct (5) 

Proof. This Q is a transport equation and we have 

i = —S(z,t), z(0) = zq 

for the characteristics z(z n ,t). As S(z,t) is odd, 5(0, t) = 0. Therefore, the origin 
is a stationary point. The estimate @ yields 

\S(z,t)\ < C\z\, Vt>0 and Vz e C 

Therefore, 

|r| <CY, r=\z\ 2 

and 

r(Q)e- ct < r{t) 

This implies ([5]) as the patch T(t) is convected by the flow. □ 
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On the PDE level, we will construct an approximate solution to the 2d Euler dy- 
namics of contours where the error can be interpreted as the strain. It is important 
to mention here that two centrally symmetric patches can not approach the origin 
in the rate faster than double exponential even if one places them in the strain S 
from the theorem ll.il In other words, the following estimate holds true 

dist(O(t),0) >exp(-e 9t ) (6) 

with some positive constant q which depends on the initial configuration. This 
is an immediate corollary of the lemma 8.1 from 2 , page 315. The constant Si 
from the theorem as well as q in ^ can be changed by a simple rescaling (i.e., by 
multiplying the value of vorticity by a constant) thus the size of Si is small only 
when compared to the parameters of the problem. 

In M 2 , in contrast to T 2 , the kernel of A -1 is easier to write and A -1 can be 
defined on compactly supported L 1 functions so we will address the problem on the 
whole plane rather than on T 2 . On the 2d torus, similar results hold. 

The interaction of two vortices was extensively studied in the physics literature 
(see, e.g., [El HI])- For example, the merging mechanism was discussed in [10] where 
some justifications (both numerical and analytical) were given. In our paper, we 
provide rigorous analysis of that process and obtain the sharp bounds. 

In [7], the authors study an interesting question of the "sharp front" formation. 
Loosely speaking, the sharp front forms if, for example, two level sets of vorticity, 
each represented by a smooth time-dependent curve, converge to a fixed smooth 
arc as t — > oo. Let the "thickness" of the front be denoted by Sf ront (t). In [7J, the 
following estimate for 2d Euler dynamics is given (see theorem 3, p. 4312) 

6front(t) > e^ M+B ^ 

with constants A and B depend only on the geometry of the front. The scenario 
considered in our paper is different as the singularity forms at a point. 

The idea of the proof comes from the following very natural question. Consider 
an active-scalar dynamics 

6 = V0- V x {Ad) (7) 
where A is the convolution with a kernel Ka(Q- I n some interesting cases, A = A~" 
with a > so Ka is positive and smooth near the origin. It also obeys some 
symmetries inherited from its symbol on the Fourier side, e.g., is radially symmetric. 
Perhaps, the most interesting cases are a = 1 (2d Euler), which is treated in this 
paper, and a = 1/2 (SQG), which is somewhat less established in the physics 
literature but still is very interesting mathematically. If one considers the problem 
with A = A~ a on the 2d torus T 2 = [— 7r, tt] 2 , then there is a stationary singular 
weak solution (the author learned about this solution from [5]), a "cross" 

O s =Xe + Xe- - XJ ~ Xj- (8) 

where E — [0,7r] 2 and J = [— 7T } 0] x [0,tt]. One can think about two patches 
touching each other at the origin and each forming the right angle. This picture 
is also centrally symmetric. Now, the question is: is this configuration stable? In 
other words, can we perturb these patches a little so that they will converge to 
the stationary solution at least around the origin? The flow generated by 6 S is 
hyperbolic and so is unstable. However, a suitably chosen curve placed into the 
stationary hyperbolic flow will converge to a sharp corner as time evolves. The 
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problem of course is that the actual flow is induced by the patch itself and so it will 
be changing in time. That suggests that one has to be very careful with the choice 
of the initial patch to guarantee that this process is self-sustaining. Nevertheless, 
that seems possible and thus the mechanism of singularity formation through the 
hyperbolic flow can probably be justified. We do it here by neglecting the smaller 
order terms. In general, the application of some sort of fixed point argument seems 
to be needed. Either way, this scenario is a zero probability event (at least the way 
the proof goes) if the "random" initial condition is chosen. However, if one wants 
to see merging for a long but fixed time, then this can be achieved for an open set 
of initial data so from that perspective our construction is realistic. 

We will handle the case a — 1 only without trying to make the strain smooth. 
The issue of whether one can choose S(z,t) 6 C°° or even drop it completely will 
be addressed in the separate publication and will, perhaps, require the realization 
of the full fixed-point argument idea explained above. We were able to do that so 
far only for the model equation where the convolution kernel is a smooth bump. 

For a > I, the same argument shows that the merging happens with exponential 
attraction to the origin. This case is much easier as the convolution kernel is not 
singular anymore so, for example, one can take the strain exponentially decaying 
in time. 

For a < I, we expect our technique to show that the contours can touch each 
other in finite time thus proving the blow-up. This, however, will require serious 
refinement of the method. 

There are other stationary singular weak solutions known for 2d Euler dynamics 
and for other problems in fluid mechanics. It would be interesting to perform 
analogous stability analysis for each of them with the goal of, e.g., solving the finite 
time blow up problem. 

The structure of the paper is as follows. In the sections 2,3, and 4, we prove 
some auxiliary results. The section 5 contains the construction of £l(t), T(t) and the 
proof of the theorem II. 1 1 In appendix, we find an approximate self-similar solution 
to the local equation of the curve's evolution. 

2. Velocity field generated by the limiting configuration 

In this section, we consider the velocity field generated by one particular pair 
of patches which resembles the limiting (t = oo) configuration around the origin 
of the dynamics described in theorem If .11 In R 2 , take 9 s (z) = xe{z) + Xe-( z )> 
where E = [0, 1] 2 . Recall that in the case of T 2 analogous configuration ([8|) is a 
steady state. Let x,y > 0, the other cases can be treated using the symmetry of 
the problem. We have 




Integrating, we have for the first component 
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where r , 1 ( 2 )(^,y) are odd and smooth around the origin. Integrating by parts, we 
get the following expression for the integral 

y log(a: 2 + y 2 ) — 2y + 2x arctan(y/x) 

By symmetry, for the second component of the gradient we have 

x\og(x 2 + y 2 ) - 2x + 2yarctan(a;/y) 

Thus, the velocity is 

u s (z) = V ± A- 1 6 S = C(-x, y) \og{x 2 + y 2 ) + r{x, y) (9) 

around the origin, where \r(z)\ < \z\. By scaling the time in Euler dynamics, we 
can always assume that C = 1/2 in ^ so we get 

Lemma 2.1. Around the origin, we have 

u s (z) = W ± A s (z) + r{z) (10) 

where 

A s (z)=xylogQ(x,y), Q(«. !/) = { |^ $ > j*j (H) 

and 

\r(z)\ < \z\ (12) 

Proof. For x,y > 0, this is a direct corollary of the calculations above. For the 
other cases, this follows from the symmetry of configuration. □ 

Remark 1. The function V^Ag is not continuous on the diagonal. That, 
however, is not so important as we are trying to control not smoothness of the 
patches but rather their distance to the origin. This particular choice of A s was 
made only to simplify the calculations below. 

We will focus later on the first term in (flO|) as the second one has a smaller size 
and will be absorbed into the strain later (see the formulation of the main theorem). 
The level sets of A s around the origin are hyperbolas asymptotically. Indeed, take 
x = ex, y = ey and consider 

A s (z) = e 2 loge 

or 

^/ + logQ(x ) y) \ =1 (13) 



loge 



For the velocity, 



W -L A /„\_J (-x(logy+l),ylogy), y>x>0 . , 

V (-x\ogx,y(logx + l)), 0<y<x 

For other x and y, we can use the symmetry of A s . 

These formulas show that the flow generated by A s is hyperbolic around the 
origin. Moreover, the attraction and repelling is double exponential. For example, 
the point (0,yo) will have a trajectory (0,yg ). In the next section, we will study 
the Cauchy problem associated to this flow. 
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3. The Cauchy problem for the dynamics generated by the limiting 

configuration 

The construction of the patches and the strain in the main theorem will be based 
on the calculations done in this section. We consider the following Cauchy problem 
where t is a parameter: 

i(r,f) = V x A s (2(r,*)), r>0 

the formula for the right hand side is given in (|14[) . and the initial position is 

z(0, t) = (e(f), tT 1 ), e(f) = exp(-e 5t ) (15) 

For now, S is some fixed positive number. So, at each time f, a point (let us say 
it has an index t) with initial position (e(i),e _1 ) starts moving under the flow so 
that at time t + t (i.e., time r spent since the start of the motion) we will see the 
arc built by points that are indexed by f + t\,ti G [0, r]. Let us call this arc Ti(t) 
(rotated by 7r/4 in the anticlockwise direction, it will be a part of the T(t) — dil(t) 
from the main theorem). 

Our goal in this section is to study the Ti(i) around the origin when t — > +oo. 
This is a straightforward and rather tedious calculation which we have to perform. 

For each t, consider the trajectory z(t, t) within the fixed neighborhood of the 
origin: {[0, e^ 1 ] x [0, e -1 ]}. This trajectory will be an invariant set for potential A s 
that corresponds to A s = —e~ 1 e(t). The parameter t, again, can be regarded as 
the tag of this trajectory. For t — > oo, around zero these trajectories will look more 
and more like hyperbolas due to (|T3l) . 

There will be two important events for each trajectory: the first is when it crosses 
the line y = x at time r = T\{t) and the second is when it crosses the line x = e _1 
and thus leaves the domain of interest < x < e _1 ,0 < y < e _1 . We will denote 
this time by 22(f). Since the trajectory is an invariant set for A s , 

z(f,T 1 (f)) = (e(f),e(f)) 

and e(t ) can be found from the equation 

? 2 (f)log?(t) = -e- 1 e(i) (16) 

This gives 

m = ^^(l + m, t-oo (17) 

One can write an asymptotical expansion up to any order but we are not going to 
need it. 

Within y > x > 0, we have y = ylogy, y(0,t) — e _1 and so 

y(r,f)=e- e \ r<Ti(f) 
Then, the value of T\ it) can be found from 

e =e(f) 

or (due to (^6| ) 

2e Tl - Ti = 1 + e st 

That gives us an asymptotics 

2\(f) = St -\og2 + e- St (l -\og2 + St) + 0(t 2 ) e - 2St (18) 
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Therefore, at the actual time w, 

w =t + T x (t) = t(l + 8) -log2 + e- 5 *(l - log 2 + St) + 0(t 2 ) e - 2St 

the curve we study will intersect the line y = x at the point (e(t),e{t)). Then, 
assuming that t is a function of w, we have 

t(w) -JL_ + Ml _ 2-^(^ + 1-106 2 + ^106 2) JltB , H 
1 >~l + S + l + S 1 + S e +U(e ) 

S 1 = S/(l + 5)<l (19) 

and 



f{w) = e(t(w)) = F{S)exp f- e ^ w 2- 1 ^ 5+ ^ + M^1_I2 W ^ (l + o(l)) 

/log 2 1 <Si . , n <S?log2 
^(*) = exp -|- - - + - - ft log 2 1 



(20) 



2 2 2 ° 2 

3+. The asymptotical form of the curve within e^ 1 > y > x > 0. We first 

compute the scaling limit of the curve in the ||z| < (f{w)j j neighborhood of 

zero where k is any small positive number. Fix the time w and solve the system 
backward in time using auxiliary functions a(r, w),0(t, w): 

d = a(log/3 + l), /3 = -£log/3 

and a(0) = f(w),j3(0) = f(w). We get then 

p{r,w) = {f{w)Y\ a(T,t) = e T (f(w)) 2 ~ e - T (21) 

Take T, the moment at which the curve is studied, and let 

T = w-r, a = f(T)&, = f(T)0 

Given a and /3, the parameters w and r are now the functions of T, S, and j3. 
We have equations 

e T {f{w)?-°- T = f(w - r)a, (/M) 8_T = f(w - T )p 

Denoting u = log a,v = log /3, we get 

_ T log/O-r) 

e 

u _ log/(w) 



w o _ T , T log/(w - r) 

z — e 



log/H iog/H 

Now, let r € [0, C] be fixed and take u> — > +oo which is the same then as taking 
T — > oo. From (|20|) . we immediately get 



v e~ T - e~ 5lT 



~ =X7 = -s-(T,ft), w^O 

u 2 — e T — e ° lT 



First, notice that 

s_(r,ft) >0 

and 

1-ft 



S _(r,ft)^o;, w= — — ie(0,l), r^O 
1 + ft 
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Now, we see that, depending on r > 0, we have different zones for x and different 
asymptotical regimes. For example, if r — > 0, the limiting curve is 

y = x~ u (22) 

Remark 1. For the distance to the origin we have 

minfx 2 + x~ 2uJ \ =xL+ x m 2ul < 2 

x>0 \ 1 

where 



Notice now that 8\ — > as 6 — > and 

w -+ 1 if 5^0 

so in the 6 — > limit the limiting shape for r <C 1 is hyperbola, which is approx- 
imately the invariant set for A s . This makes the perfect sense as the initial point 
(e _e , e _1 ) approaches the separatrix OY slower when 6 — > 0. Had it been constant 
in t, the curve would exactly follow the invariant set for A s . 

Where precisely on the rescaled space x, y will this uniform convergence take 
place? We have 

a(r, t) = e T (/(T + r)j ~ (/(T)j « f(T) 

where Ti(T) — > arbitrarily slowly. 

For each r > 0, there will be its own part of the curve with the corresponding 
scaling limit given by 

y = r s -W'' 

For this fixed positive r we can easily identify the domain where this scaling is 
going to take place. Indeed, at time T, the formula (|2"0)) yields 



(23) 



f(T) - F(5)exp ^-e^-W+i) + S -±^L ^ rj (1+5(1)) 

For the j/-coordinate of the point on the curve in the domain represented by r, (|2T 
gives 



y 



exp | 



_ 2 -l/(l+<5) e <5iT 



di(r) 



exp ( ^ T 



where di(r) = e"( 1 " ,5l)T < 1. Thus, we see that y - (/(T)) dl < T ) ± and for the 
^-coordinate we have similarly 

x - (f(T)) h(T)± , h(r) = (2 - e~V lT > 1 

For any fixed r > 0, we computed the asymptotical shape of the curve. However, 
we will also need the bounds on the curve for t(T) — > +oo. Given arbitrarily small 

fixed fi > 0, the the part of the curve in the region y > f/(T)^ will be contained 

in the rectangle 

0<x<f(T), (/(T))" <y<e- 1 (24) 

The first inequality is obvious and can certainly be improved a lot but this is all 
we will need later on. 
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3.2. Asymptotical form of the curve within < y < x < e 1 . This part of 
the curve is more complicated. To study it, consider equations 

d(r) = -a(r)loga(r), 0(t) = /3(r)(log a(r) + 1) 

with initial conditions a(0) = /(w),/3(0) = f(w). We are interested in the shape 
of the curve at time T — w + r. 

a{r) = {f{w)y-\ = (/H) 2 " 6 " V (25) 

We again scale by /(T) as follows 

a = f(T)a, f3 = f(T)p 

If u = log a, v = log (3 then 

T +2 _ e -T _ l °Sf( W + T ) 



v = log f(w) log/(w) 

u e _ T _ log/(w + r) 

log /(to) 

7) 9 — P~ T — P SlT 

We again have 

s+(r,5i) -> w 

so the rescaled curve converges to the graph of 

/ \l-ri(T) 

and this convergence is uniform on f(T) < x < (/(T)l where n(T) is 

positive and tends to zero arbitrary slowly 
For any fixed re (0, r cr ), we have 

2- e- T - e 5lT 
e~ T - e 5 

The critical T cr {8\) = log^ cr {8\) where £ C r(#i) is the solution of the equation 



y = x 



*+Mi) = — — oSlT >o 



2-1-^ = 0, e > i 

Notice that £ cr (<5i) —¥ 1 as <5i — > 1 and £ cr (#i) -4- +oo as 8\ — > 0. Consequently, for 
the critical value 

TcrOl) -> 0, 5l -> 1 

and 

T cr {Sl) -> +00, #1 -> 

As before, for the subcritical value r < r cr , the curve will have a limiting shape 
given by 

however s+(r, <5i) decays in t and s + (r cr , <5i) = 0. 

The (a;, y)-coordinates of part of the curve corresponding to each r < r cr is easy 
to find. We again have 

f(T) = F(S) cxp ^- e ^ T 2- 1 /(^+i) + W ~ 1) T ^ (1 + 5(1)) 
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and 

x ~ (/(T))^)* d 2 (r) = e -^ 1+5 i) < 1 

ia(T)± = 0f ,-5iT _ e - T (l+«i) 



y ~ (/CO) > « 2 (r)=2e 



and ^(t) > 1 for r < r cr . In the x-coordinate, the domain corresponding to r < r cr 
will be characterized by /(T) < x < f d ^( T ^-)(T). Moreover, 

What can be said about the curve for the region with r > r cr ? For the analysis 
that follows in the next section we will only need very rough bounds. 

Starting at time w, the point with coordinates (f(w),f(w)) will go along the 
trajectory which will then cross the line x — e _1 in time 

t = log log —3— - Siw- +5(1) 
f(w) 1 + 6 

From we get for r > r cr and large w 

R f \ 2-2e~ T / N 2-2e~ T <='' / \ 2-2e _Tcr 

£ = ^(/Mj < (/Hj < (/(T/(l + <*!))) 

Therefore, part of the curve corresponding to r > r cr will be contained in the 
rectangle 

[(f(T)) MTcr ~\e- 1 } x [0, (/(r/(l + <5 1 )) <(5l) ], aSi) > (27) 

3.3. Self-similar behavior around the origin. Summarizing the calculations 
above, we conclude that around the origin the rescaled curve will converge to the 
graph of the function y — . How large is the domain of convergence? If the 
curve is calculated at time T and if T\ (T) is arbitrary positive function such that 
T\(T) — ¥ as T — ¥ co, then this convergence will be valid in the region 

f(T)) <x<(f(T) j 



Outside this window, we have different scaling limits for different values of the 
parameter r > 0. 

Remark 2. What is the nature of the scaling law for r ~ that we have got? 
Well, the Euler dynamics is defined by the convolution with the kernel 

K(z)~]ag\z\ 

and | log(se)| = | loge| + 0(1) for s G [1,M] with any fixed M. Thus, the strength 
of the created hyperbolic flow is more or less the same within any annulus e < \z\ < 
Me. Outside this annulus, say, for \z\ = yfk, the size of the kernel is quite different. 

Now, assume that we are given the hyperbolic flow in the whole plane defined 
by the following equations 

x = y 

y = x 

Let us find the evolution of the contour C(t) under this flow which would be self- 
similar in the sense that 6(f) = Q(0)e~ 5t . Assume that C(i) is given by the graph of 
the function y(x, t) and then one gets nonhomogeneous Burgers equation for y(x, t) 

y t (x, t) = -y x (x, t)y(x, t) + x 

By our assumptions, 

y(x,t) = e- St H(e St x) 
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so we have 



(28) 



H' = 



5H + Z 
5£ + H 



IfH = £z(£), then 



\z-l\ A \z + l\ 



B 



c\t\ 



1 + 5 



2 



,B = 



1-6 
2 



and so 



£,\ A \H + i\ B = C 



Due to scaling, we can assume C = 1. In the coordinates H — £ = /3, H + £ = a, 
we have 



(compare with (|2"2"j) 1. 

Is it possible to find the initial data such that the evolution of the curve is self- 
similar on the larger interval? The answer to this question is yes. For example, in 
(p~5|) . one can take 



with 7 € (0, 1). Then, the distance from the curve to the origin will be < expf — e* 



71 G (0, 1) but the self-similar behavior will take place in a wider relative range (but 
yet not on the whole ball B e -i(Q)\). This is a general rule: the slower the Cauchy 
data jumps from one invariant set of A s to another, the more regular the curve is 
around the origin. 

4. Comparison of the velocity fields generated by T 1 (t) and by the 



Let us recall that we denoted by T± (t) the arc constructed in the previous section. 
Let us take Fi(t) and close it in the smooth and arbitrary way to produce the 
simply-connected domain Cli(t) (see Figure 1). We will compare now the velocity 
fields generated by two patches E U E~ and fii(i) U fij"(f). Recall that 



j3 = a 



(l-«5)/(l+5) 





LIMITING CONFIGURATION. 





O 



Figure 1 
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We will need to use the following lemma in which the picture above will be rotated 
by 7r/4 degrees in the anticlockwise direction. Let us denote by ^(t), ^{t), and 
Ei the resulting sets. 

Fix t and assume that some arc (let us call it T 2 as we will later apply this 
lemma to the part of ^(i)) lies above the graph of |x| and below the graph of a 
certain function g(x) defined on [—a, a] where a = l/(ev / 2)- In this section, we take 
e = /(*)■ 

Lemma 4.1. Rescale g(x) as follows 

g(x) = e~ 1 g{ex), \x\ < a/e 
Assume that g satisfies the following properties: 

(a) < C\ < g(x) < C*2, x G [—1, 1] with some absolute constants Ci( 2 ). 

(b) For \x\ > 1, 



\x\ < g < \x\ + h(x), h(x) = o(\x\), x — > oo 

Consider 

D{z,t) =W ± A- 1 ( X(M t)un-(t) ^ Xe 2 ue-) 
Then, we have the following estimate for \z\ > e 

\ 



(29) 



|I>CM)I<I*I 



\ K|x|<o/e 



h(x) + h(-x) 
\x\ ■ (\x - zA + 1) 



-dx 



z = ez,z = (zi,z 2 ) 



) 



Proof. Take z G C + . From the central symmetry, we have 
V ± A- 1 ( x ^ uO --x B2UiS -)(0)=0 
Now, let us use the following formula 



\i-A 2 \H? 



to get 



\D(z,t)\<\z\ 



1 



\ (B 2 \n 2 )n{|z|<a/io} 



l*-£H£l 



+ 



+ 



/ 



(E 2 \n 2 )-n{\ z \<a/w} 



for any e < \z\ < a/10. Scale by e and notice that \z — £\ > \z\ — £i| so the first 
integral can be estimated by 



C 1 



h{x) 



-dx 



H<\x\<a/e \x\ ■ (\X-Z!\ +1) J 

For the second integral, we have a similar bound with h(—x) by symmetry. 



□ 



We also immediately get 
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Corollary 4.1. If the function h in i29\) satisfies 

\h(x)\<\x\-*, 7 >0 

for \x\ < e~ Q , a e (0, 1) and 

h(x) 



| In e| max 

e _a <|5T|<a/£ 



< c 



(30) 
(31) 



then 



Proof. From O 



and (jnU) yields 



\D(z,t)\<\z\, \z\>e 
h(x) 



K|x|<£-« 1^1 ' + !) 



(lzil + 1)- 1 



— cJx < line] max 

:\x\<a/e Fl ' IF - 21 1 + !) £-"<|£|<a/e 



h(x) 



< C 



□ 



Now, let us apply this corollary to our situation. For this, we will use the results 
from the previous section. 

Remark 1. Notice that for each 6, the arc Ti(t) satisfies the conditions of the 
corollary with some a(5). For |a;| > e~ Q , we can use (f24|) and (|27|) to get 



h(x) 



< (/(V(1 + ,5 1 ))) C , C>0 



and 



Thus, we have 



(log/(t))/(t/(l + *i)) 



sup 

M>/(*) 



D(z,t) 



l + 5i 



< C 



(32) 



Remark 2. For \z\ < f(t), we can argue as follows. Instead of E above, consider 

E(t)=E\{\z\<f(t)} 

Then, similarly to (|32]l . 

V^A 1 (Xn 2 (t)un-(t) _ Xe 2 (*)□£-(*)) ~ M 

but now this is true for all z S £?i(0). For V" L A _1 XB 2 (t) , the direct calculation 
yields 



V^A~ 



l X E2{t)uE - (t) =CV x (xylogf(t)) +0(\z\), \z\<f(t) 



and so 



2 (t)un 2 -(t) 



CV 



- (ay log 0(H), |*|</(t) (33) 



where we again can assume that (7 = 1. 
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5. Construction of the vortex patch dynamics and proof of the main 

THEOREM. 

We first construct an incompressible strain ^Sf which satisfies the following prop- 
erties (see Figure 2 for the upper part of the actual picture): 

1. ^ is odd and is compactly supported. 

2. Around the points (4,4) and (—2,2) it is the standard hyperbolic time- 
independent flow (these are the domain D\ and Dq). At (—2,2), we choose the 
separatrices to be y\ = — x and j/2 = 2 + (x + 2) . The flow is attracting along ?/2 and 
is repelling along y\. To define the curve T(t) in D\ we need the following result 
which will later guarantee necessary initial conditions for the dynamics around 
(0, 0). Let the local coordinates near (—2, 2) be denoted by (£, n). 

Lemma 5.1. Fix any 8 > and consider the standard hyperbolic dynamics around 
the origin 

i = i, c(o) = & 

V = -V, 17(0) = Vo 

Let G(£) = for £ < and G(£) = ^ exp(-^ s ) for £ > 0. Consider the 
evolution of the smooth curve T(0) = {(£, G(£)), |£| < 1} under this flow. Call it 
r(i) = {(£,G(£,t)),|e|<l}. Tften, 

G(l,t)=e- e " (34) 

Proof. This is a straightforward calculation. The point (£o, Vo) moves to (£oe > ?7oe _ *) 
in time Thus, the point (e _t , G(e _t )) G T(0) will move to (1, e~ e ) in time t. □ 

Remark 1. Notice that the part of T(t) that belongs to the left half-plane does 
not change in time. Within the window |£| < 1, the curve T(t) is always smooth 
and converges to the coordinate axis. 

3. Around the origin (domain D2), we choose 

* (2 , i) = V -A,(^,^, t ) (35) 

where 

A s (z), z£l-0.01f(t) 7 0.01f(t)} 2 
xylogf(t), ze [-0.01/(i),0.01/(i)] 2 

Recall that A s is defined in (fTTj) and f(t) is introduced in ([23)1 . We changed the 
coordinates in ([35]) as we want to rotate the picture described in section 3 by 7r/4 in 
the positive direction. We also modified the value of the potential in the 0, 01/(i) 
neighborhood of the origin to get rid of the artificial singularity generated by the 
sharp corner in the limiting configuration. Notice now that 

*(z,t)=(^ + { l))^(yz-x*) 

in {\z\ < mf(t)} with any fixed m. 

4. Between D\ and D2 the potential can be smoothly interpolated. 

5. In ^3(5) , the flow is laminar with direction perpendicular to the black segments 
and in the north-eastern direction. 



A x (z,t) 
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6. The potential between zones D2 and D3 can be smoothly interpolated, as 
well as the potential between D5 and Dq. In the zone Dj, the potential is zero so 
the curve is frozen. This zone again is interpolated smoothly between D\ and D$. 

7. In the zone D4, we construct non-stationary potential in the following way 
(only in this zone the flow is essentially time-dependent!). We need an argument 
that allows an interpolation between two laminar flows and guarantees the pre- 
scribed evolution of the curve in these laminar zones. What we want is to 
define dynamics in the regions D3, D4, D§ right after the points on the curve leave 
Z?2- We need to define this dynamics in such a way that the motion of T(t) is 
localized to these regions and, moreover, that it does not move in D5. Once again, 
in D3 and -D5 we postulate the flow to be laminar and then we want to define it in 
D4. We will do that in the local coordinates. 

Assume that potential A(z) = — y in B = {z : — 1 < x < 0} U {z : 1 < x < 2} ~ 
D3 U D5 . This potential generates the laminar flow 

6 = Vtf-V^A 

where V^A(z) = (—1,0). We want to define smooth A(z, t) in D4 — {z : < x < 1} 
such that the resulting A(z, t) is smooth globally on D3 U D4 U D5. Moreover, given 
smooth decaying v(t) (e.g., v e i 1 (R + ) is enough for decay condition), we need to 
define a curve T(0) = {(x, -f(x, 0))} that evolves under this flow T(t) — {(x, 7(2;, t))} 
such that 7(0,4) = v(t) and 7(1, t) = 0. This function v{t) is determined by T(t) 
in the zone D2 where it approaches the separatrix in the double exponential rate. 
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To be more precise, v is proportional to the distance from T(t) to this separatrix 
in the area where Z? 2 and D3 meet. 
We will look for 

AO, t) = -y - gi (x)g 2 (x - t) 

where 51(2) are smooth. Then, to guarantee the global smoothness, we need gi(x) — 
around x — and x = 1 . Now, take a point (0, v(T)) and trace its trajectory for 
t > T. We have 

x(t,T) =t-T,te[T,T+l] 

and 

y(t,T) = v{T)-J (<?i(T-r) 52 (T-T-T)+.g 1 (T-r).9 2 (T-T 1 -r))dT, * e [T.T+l] 
Since we want y(T + 1,T) = and g\ to vanish on the boundary, 

v{T)=g' 2 {-T) f + gi (T-T)dr = g' 2 (-T) f 9l (x)dx 

JT JO 

and this identity should hold for all T > 0. Take any gi with mean one, this defines 
52 on the negative half-line as long as we set #2 (—00) = 0. We can continue it 
now to the whole line in a smooth fashion to have g 2 globally defined. How do we 
define the initial curve at t = 0? We extend smooth v{t) to t £ [—1,0] arbitrarily 
and apply the procedure explained above to t G [— l,oo). The curve that we see 
at t = will be the needed initial value for the dynamics that starts at t = 0. It 
is only left to mention that to localize the picture in the vertical direction we can 
multiply A(z, t) be a suitable cut-off in the y direction. 

The part of the curve that is in D5, Dq 7 Dj, and the north-western part of D\ is 
stationary, it does not move at all (this is easy to ensure by making this part of the 
curve the level set of the stationary potential A(z), \& = V^A). For the rest of the 
curve, it does change in time and the flow is directed along it in the anti-clockwise 
direction. 

Now that the explicit ^(z,t) and the curve T(t) = dft(t) evolving under this 
flow are defined, we are ready to prove theorem 1 1.11 

Proof. (Theorem \l.l\) . We have by construction 

9 = V6- $>(z,t) 
and 9(z,t) = Xfi(t)O) + Xn-(t){z)- Let us define S(z,t) by 

S(z,t) = ^(z,t)-V ± A- 1 8 

Then l[52j). and jnj) give © for S(z,t). (Here we also observed that the 

contribution from UO _ (t))\5 .i o (0) to V^A - ^ is of order 0(\z\) around the 
origin^ . 

Thus, around the origin, the error S{z,t) will in fact be much smaller than 
i$f(z,t). It is odd as is even and is also divergence free as the difference of two 
divergence free vector fields. We get 

= ve • (y^A^e + s(z, t)) (36) 
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and the theorem is proved as the dynamics of f2(t) satisfies 

dist(fi(i),fr(i)) = 2dist(n(i),0) ~ f(t) 
For f(t) we have a double exponential decay due to ([23]) . 

□ 

Remark 2. Strictly speaking, the Cauchy data for the evolution in the domain 
Z?2 will not be given by ()34|) as the flow will distort it when moving between zones 
Di and D^. However, this leads only to a minor change (a fixed time increment, in 
fact), so the same argument goes through. 

Remark 3. The vector-field ^ and the strain S we constructed are not even 
continuous, e.g., ^ jumps on the lines y = ±x. That, however, is not an issue 
as the dynamics can still be defined and we are controlling the size of S, not its 
smoothness. In fact, we could have chosen A s = xylog(x 2 + y 2 ) in (JXTJ) . That takes 
care of the regularity issue but the calculations in section 3 would have been much 
longer. 

Remark 4. The parameter 8\ in the formulation of the theorem ll . ll corresponds 
to 8\ in (fT9| and can be chosen arbitrarily from the interval (0, 1). The size of this 
interval is determined by the parameters of the problem: the value of vorticity, 
initial size of the patches and their distance to the origin. 

6. Appendix: approximate self-similar solution to the contour 

dynamics 

In this section, we address the following question: is it possible to construct the 
approximate solution to the Euler dynamics of patches such that the self-similarity 
will persist on the larger set? We will do that in a rather artificial way as we already 
can make the right guess for what the solution should be. 

Assume that the boundary of the simply-connected patch is parameterized by 
j(s,t). Then, the velocity at every point of the contour can be computed by ([5], 
formula (1)) 

p2TT 

u(7(£ t),t) = C log | 7 (f , t) - 7 (s, t)\j' s ( s , t)ds 
Jo 

This is a simple corollary of the Gauss integration formula. 

If we have a centrally symmetric pair of vortices interacting with each other and, 
like before, the part of T(t) close to the origin can be parameterized by the function 
y(x,t), then the equation for evolution reads 



(37) 

where r(x, y, t) is a contribution from those parts of T(t) and T (t) that are away 
from the origin. This r(x,y,t) is therefore smooth and r(0,0,t) = by symmetry. 
Let us drop r(x, y, t) and try to find an approximate self-similar solution? In other 
words, we want y(x,t) — e(t)</>(x/e(t)) to satisfy (|3T[) up to some smaller order 
correction. Substitution into (|37p gives 

- (4>(x) - cf>'(x)x) = (38) 
e 
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0.5<T ' 



J \ (x + £ ) 2 + U x ) + d> £ 2 / 



— 0.5e _I 

where E is an error we will control later on. Let us rewrite the integral as follows 
(up to a constant multiple) 

{<j>\x)-^\mn+<i>{x)m)K{x,M 



-0.5c- 1 

with 

a = (x - £) 2 + (J>(x) - <H£)) 2 , b = (x + £) 2 + {4>{x) + 
and so we get 

k 1 (x,t)(p\x)x + k 2 (x,t)<j>'(x)<j>(x) + k 3 (x,t)x + ki{x,t)<j){x) (39) 

where the coefficients kj are defined correspondingly. Assume now that <f) satisfies 
the following assumptions: 

(a) <f)(x) is smooth 

(b) < Ci < (f>(x) < C 2 for x € [-1, 1] 

(c) 4>(x) = \x\ + pix) where |x| > 1 with 

\p(x)\<\x\-\ \p'{x)\<\x\-\ 7 >0 

Let us estimate the coefficients kj now. We will handle k 3 , the analysis for k 2 is 
similar. 

Consider x G [f , 0.56" 1 ], the other values can be treated similarly. For the integral 
over the positive £ we have (after the change of variables £ = x£i , recall that x = ex) 

5x — 1 

J a 1 + (i+p(x)/mi+pmi)/x) ^ 





I /-(i+ei) 2 + (i+«i+p(£)/£+p(£?i)/£) 2 drj 

4 J(ii~i) 2 +(i-ii+ P (sys- P (si;i)/s) 2 V 



For A, we have a representation 



so 

0.5x~ 



1 

= 0.51oga; + O(l) 



THE SHARP CORNER FORMATION IN 2D EULER DYNAMICS OF PATCHES 



19 



For the other integral, changing the sign in integration 

°- 5e 1 b 

h = - I £(l + p'(-0)r^- / 

J b-a J a T) 



o 

Here, we have 

b = (x - 2 + (x + £ + p(x) + p(-0) 2 , a = (x + 2 + (x-t + P(x) - p(-0) 2 
As both x and £ are large in the interesting regime, we are in the situation when 

a 7 b> (x 2 +£ 2 )/2 
so we can use the mean- value formula 

1 f b dn 1 b — a . . 

b-aj a rj b 277f 
Substituting, we have two terms: I2 = — (Ti + I2). 

<-0.5a: 



T 1= / ^(i+p'H^i))^- 1 ^! 
Jo 



where 

B = 2(1 + £ 2 ) + 2(1 + ii){p{x)/x + P (-x^)/x) 
+ (p(x)/x + p(-x^)/x) 2 

Thus, 



~ 6( 1 + ^(-^ 1)) 
2(i + d 2 ) 



T 1= / ^7,?*"" ^ +0(1) = -0.51ogs + O(l) 



For the other term, we have 

~\ x\p(-m + 

(x^+e) 2 



rp , < , x\p(-0\+Z\p(x)\ + \p(x)p(-0\ 

Jo 



S r^(i + x + ^< c 

Combining all terms, we have 

k 3 = logx + 0(1), x > e 

For a; <~ 0, we get I^) = 0.5 log e + 0(1). These calculations show that 

_ f log N+0(1), N>e 
K3 ~\ lege + 0(1), |z|<e j 

Analogous estimates can be obtained for fc 2 . They yield 

_ / log|x|+0(l), \x\ >e 
2 " I loge + 0(l), \x\<e 

The estimates for other terms are 

1*1(4)1 =0(1) 

Indeed, 

fci=/°' 5e e^oc fc 4 --/ ' 5e 0'(o^(o«"(^Ode 
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and if one does the same analysis as we did for k^ in (|40[) . we will get the sum of 
two integrals: one over positive £ and the other one over negative £. Each will have 
the same large logarithmic leading term but they will come with different signs now 
and so will cancel each other in the sum leaving us with the uniformly bounded 
error terms. 

Thus can be written as 

-loge(»-£) + E 

where 

E = (k 2 + log e)(f)'(j) + (&3 — log e)x + k\X(j)' '(x) + k±cf)(x) 

Going back to (|38|) and choosing C appropriately (C < 0), one wants to make the 
following choice for e and <f>: 

e 1 = eSi log e, x - <j/<f> = -5i(<f> - x<f> f ), Si € (0, 1) 
Take for e(t) one particular solution 

e(t) = cxp(-e' 51 *) 
The equation for <f) we had before (see (|28|) ) and so we have 

- x\ A \<j) + x\ B = 1, A=(l + <5 1 )/2, B = {l-5 l )/2 

and its solution trivially satisfies assumptions (a), (b), (c) mentioned above. 

For the original equation (|57|) . the error one gets after substituting e(f>(x/e) 
amounts to eE where 

eE < e if \x\ < e 

and 

\eE\ < e + e\(log\x\ -loge)(^-x)| + \x<f/{x)ki\ + \e<j>(x)k A \ if |ar| > e 
Therefore, for > e, 





X 


-7(<5i) 


X 


z\+e 




log 






e 




e 



and z — (x, e<p(x)). Thus, we see that the error is small again so it is possible to find 
the approximate solution with the self-similar scaling that holds on the ball of size 
~ 1. We are not trying to make this picture global and define the incompressible 
strain on the whole R 2 which corresponds to the error eE but we believe it is 
possible. 
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